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Introduction 

The reacion mechanism behind multi 
fragmentation in heavy-ion collisions has 
been a long standing question in intermediate 
energy nuclear physics. Not only the mul- 
tiplicities of various kinds of fragments are 
useful observables, but how and when these 
fragments are realized in phase space can 
also shed light on reaction dynamics. These 
correlations can be of importance to reveal 
the phase space details of decaying system. 
In this direction, Puri and collaborators have 
devised a faster cluster recognition algorithm 
namely simulated annealing clusterization 
algorithm (SACA) [H-i]. This algorighm is 
found particularly successful in describing 
spectator matter break-up at relativistic 
@, IH and ultra low @ incident energies. The 
phenomenon of global universality of rise and 
fall pattern in the multiplity of intermediate 
mass fragments was accuratley reproduced 
with this model It, however, remained 
to see how these fragments are correlated in 
phase space and originate from which region 
of the nucleus. In the present paper, we shall 
compare the phase space characteristics of 
the intermediate mass fragments (IMFs) pro- 
duced using minimum spanning tree (MST) 
and SACA clusterization algorithms. 

The Model 

The phase space of nucleons is recorded at 
different time steps using quantum molecular 
dynamics (QMD) transport theory Q. In 
QMD approach, the trajectories of nucleons 



follow Hamilton's equations of motion: 

Pi = -{r^H}, (1) 
r<= {Pi,K}. (2) 

Here H is the total Hamiltonian of the system 
of Ap + At nucleons. These equations are 
solved after fixed time interval At chosen to 
be very small. The phase space of nucleons is 
clusterized using MST and SACA methods. In 
SACA model, pre-clusters obtained with min- 
imum spanning tree procedure are subjected 
to a binding energy cut [3, [|| : 
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with Eund = 4.0 MeV if N f > 3, else 
Ebind = 0. In this equation, Eu n d is the 
fragment's binding energy per nucleon, Nf is 
the number of nucleons in a fragment, and 
Pj^ is the center-of-mass momentum of the 
fragment. 

Illustrative Results 

For this study we simulate a single event 
of 197 Au+ 197 Au collision at incident energy 
Eiab — 1000 AMeV and impact parameter b = 
8 fm. We have employed a soft EoS (X=200 
MeV) along with energy dependent n-n cross 
section. We display in Fig. 1 the distribution 
of nucleons bound in IMFs [5 < A < 65] using 
SACA (at 60 fm/c) and MST (at 200 fm/c) 
methods. One can see that higher yield of 
fragments is obtained using SACA method at 
forward as well as backward rapidities. MST 
approach, on the other hand, predicts only sin- 
gle cluster of mass A = 5 in the backward 
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FIG. 1: The distribution of nucleons bound in 
IMFs is displayed in the reaction plane (i.e. Z-X 
plane) using SACA (left) and MST (right) meth- 
ods. 
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FIG. 2: Nucleons distribution in the two colliding 
Au-nuclei in the Z-X plane at the time of initial 
contact. Filled circles in red indicate positions of 
nucleons emitted as intermediate mass fragments. 



hemi- sphere. This shows that SACA method 



is much faster in recognizing the proper frag- 
ment yields from the spectator matter break- 
up. 

Next, we backtrack the original location of 
these nuclei in the colliding nuclei at the start 
of nucleus- nucleus collision. Figure 2 shows 
these nucleons emanating as bound clusters 
in the two Au-nuclei. One can see that SACA 
method leads to exaggerated yield of frag- 
ments out of spectator domain. In MST 
method, however, nucleons are confined to 
overlapped zone only. This results in under- 
estimated yield of fragment multiplicities in 
ALADiN data. 

Summarizing, we addressed in this paper 
the problem of origin of fragments via phase 
space analysis of MST and SACA clusters. 
Our model predictions highlights the impor- 
tance of binding energy criterion in recogni- 
tion of most bound fragment configuration. 
This criterion is important to understand the 
spectator matter physics and accurately pre- 
dicting fragment multiplicities. 
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